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Thermal fluctuations in moderately damped Josephson junctions: Multiple escape 
and retrapping, switching- and return-current distributions and hysteresis 
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A crossover at a temperature T* in the temperature dependence of the width a of the distribution 
of switching currents of moderately damped Josephson junctions has been reported in a number 
of recent publications, with positive da/dT and IV characteristics associated with underdamped 
behaviour for lower temperatures T < T* , and negative da/dT and IV characteristics resembling 
overdamped behaviour for higher temperatures T > T* . We have investigated in detail the behaviour 
of Josephson junctions around the temperature T* by using Monte Carlo simulations including 
retrapping from the running state into the supercurrent state as given by the model of Ben- Jacob 
et al. We develop discussion of the important role of multiple escape and retrapping events in the 
moderate-damping regime, in particular considering the behaviour in the region close to T* . We 
show that the behaviour is more fully understood by considering two crossover temperatures, and 
that the shape of the distribution and (j{T) around T* , as well as at lower T < T* , are largely 
determined by the shape of the conventional thermally activated switching distribution. We show 
that the characteristic temperatures T* are not unique for a particular Josephson junction, but 
have some dependence on the ramp rate of the applied bias current. We also consider hysteresis in 
moderately damped Josephson junctions and discuss the less commonly measured distribution of 
return currents for a decreasing current ramp. We find that some hysteresis should be expected to 
persist above T* and we highlight the importance, even well below T* , of accounting properly for 
thermal fluctuations when determining the damping parameter Q. 

[Accepted for publication in PRB; ©American Physical Society 2008] 

PACS numbers: 74.40.-|-k, 74.50.-|-r 
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1. INTRODUCTION 

The Josephson junction system has been extensively 
studied both theoretically and experimentally. Theoreti- 
cally it has been considered a model system for studying 
escape from a metastable potential well. Experimentally, 
Josephson junctions have found numerous applications 
and are presently being used in several quantum bit im- 
plementations. In such experiments, an understanding of 
the influence of thermal fluctuations is crucial in devel- 
oping applications. Josephson junctions can be charac- 
terised by a damping parameter Q. The majority of the 
large body of previous work in the literature has concen- 
trated on junctions in either the underdamped {Q 3> 1) 
or overdamped {Q ^ 1) limits. In this paper, we focus 
on the intermediate "moderately damped" (Q ~ 5) limit, 
where thermal fluctuations lead to interesting physical ef- 
fects. 

For strongly underdamped Josephson junctions under 
the influence of thermal fluctuations, the IV character- 
istics are hysteretic and the dynamics of switching from 
the zero-voltage supercurrent state to the finite-voltage 
resistive phase-slip state are well described by the anal- 
ysis of Fulton and Dunkleberger-, with a distribution in 
switching currents as a result of thermal fluctuations. In 
contrast, overdamped junctions show non-hysteretic be- 



haviour, with a finite voltage on the supercurrent branch 
of the IV characteristic, associated with thermally acti- 
vated phase diffusion, and thermal fluctuations leading to 
very much smaller variations in the switching behaviour. 
Phase diffusion in junctions with hysteretic IV character- 
istics has been discussed by Kautz and Martinis^i and is 
associated with frequency-dependent damping, such that 
junctions are underdamped at low frequencies, but in the 
overdamped limit at high frequency. 

The temperature dependence of the switching current 
and the width of its distribution are experimental param- 
eters of much recent interest. Experimental evidence of a 
crossover in the temperature dependence of the switching 
current was reported first by Franz et al.^ in experiments 
on small "intrinsic" Josephson junctions (IJJs). They 
obtained IV curves characteristic of underdamped junc- 
tions below a crossover temperature and IV curves char- 
acteristic of overdamping above that temperature. More 
recent experimental papers have reported a crossover in 
the temperature dependence of the width a of the switch- 
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at a temperature T* , with positive 
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da/dT below T* and negative da/dT above T*. This 
was associated with a regime of moderate damping. The 
negative da/dT region was associated with retrapping 
of the phase following escape. The low-temperature be- 
haviour fits the expectations for underdamped junctions, 
and the high-temperature behaviour resembles previous 
observations for overdamped junctions with phase diffu- 
sion. One might simply explain the crossover from un- 
derdamped to overdamped behaviour by a temperature- 



dependent damping Q and this was indeed the suggestion 
of Franz et al. However, it was demonstrated by Kras- 
nov et alM- that such a crossover should also be expected 
even for temperature-independent Q if the junctions are 
in the moderately damped regime (Q ~ 5). Krasnov et 
al. derived an approximate quantitative formula with 
T* = T*{Q), implying that T* is a measure of the damp- 
ing. 

Several theoretical treatments of the retrapping pro- 
cess have been presented^>^>^ and the analysis of retrap- 
ping was conducted in various ways in the experimen- 
tal reports of a crossover in a{T). In the analysis of 
Ref. [ij, retrapping was assumed to be determined purely 
by energetic considerations: retrapping is certain to oc- 
cur where it is energetically expected, below a current 
Jm ~ 4/c/7rQ, where /c is the (fluctuation-free) critical 
current of the Josephson junction. For / > /m there is an 
energy cost At/jj to retrapping — in Ref. 113 . retrapping 
was neglected for / > Im- Krasnov et alt^ treated re- 
trapping above /m as a thermally activated process, with 
an energy barrier A[/j^, using the model of Ben-Jacob et 
al.^ Mannik et al}'^ used Monte Carlo simulations with 
an RCSJ model including frequency-dependent damping 
to determine the probability of thermally induced retrap- 
ping following escape. 

In this article, in order to conduct a semi-analytic 
analysis of the multiple escape and retrapping pro- 
cesses, we have adopted the model of Ben- Jacob et 
al. We also include the effects of frequency-dependent 
damping (see Section 14. ip . We develop discussion of 
the important role of multiple escape and retrapping 
events in the moderate-damping regime and present re- 
sults of Monte Carlo simulations showing the variation 
with experimental parameters of the mean and width 
of the switching-current distribution^ We consider the 
crossover between the lower-temperature conventional 
underdamped regime and the higher-temperature over- 
damped regime. Although previous studies have consid- 
ered a single crossover temperature, we show that, in 
detail, the change occurs in two stages, with a lower- 
temperature transition from underdamped behaviour to 
behaviour in the crossover regime, and a higher temper- 
ature transition from the crossover regime to the higher- 
temperature overdamped regime. We demonstrate a sig- 
nificant change in the shape of the switching current dis- 
tribution around the crossover and study this quantita- 
tively through the skewness parameter. 

The process of return from the resistive state in a hys- 
teretic junction is a much less well-studied phenomenon 
than that of escape. Here we also consider the process of 
return from the resistive state to the supercurrent state 
as the current is ramped down, and the resulting varia- 
tion in hysteresis around T* . We compare our findings 
with previous reports in the literature. 

In our Monte Carlo simulations, for a current /, the 
probability in a short time interval 5t of a transition 
between the metastable and running states is given by 
rg(/)(5i for escape from the metastable state (with Fg 



given below by Eqn. ^ or T-n(I)5t for retrapping from 
the running state (with Fp given below by Eqn. [3]). A 
bias current is ramped up (or down) at a constant rate 
in order to generate distributions of switching (or re- 
turn) currents for junctions with a number of different 
parameters. We neglect the temperature dependence of 
the critical current Iq and Q in order to emphasize ef- 
fects due to thermal fluctuations in the junctions. As a 
bias current is ramped up, a switch is counted when the 
junction spends more than half the time in the running 
state over some time period r.^'' Throughout this arti- 
cle, we use the term "escape" to describe any (possibly 
short-lived) escape from the instantaneous zero-voltage 
state, and reserve the term "switch" to describe an ex- 
perimentally measured switch to the running state. Sim- 
ilarly, when describing the behaviour as an applied cur- 
rent is ramped down from the critical current, we re- 
serve the term "retrapping" to describe a (possibly short- 
lived) change from the voltage state to the zero-voltage 
state, and use the term "return" to describe an experi- 
mentally measured change from the voltage state to the 
zero- voltage state. 



1.1. The RCSJ model - the underdamped regime 

For a resistively shunted Josephson junction in the ab- 
sence of fluctuations, escape from the supercurrent state 
to a state of finite voltage characterized by the junction 
resistance occurs when the current bias applied to the 
junction reaches the junction critical current Iq. At fi- 
nite temperatures, thermal fluctuations lead to switching 
at currents below Iq, and there arises experimentally a 
distribution in possible values of the switching current. A 
common experimental configuration is to ramp the cur- 
rent up from zero at a constant rate d//dt. In that case, 
the probability of a switch in the current range / to /+d/ 
is p{l)dl, with^ 
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where Tt? is the rate, at current /, of escape from the 
supercurrent state. 

A mechanical analog for the resistively and capaci- 
tively shunted Josephson junction (RCSJ) is that of a 
particle in a washboard potential; it is often used in dis- 
cussing the dynamics of such junctionsJ^ The height of 
the corrugations in the untilted washboard is set by the 
Josephson energy. The current bias corresponds to a tilt 
of the washboard, and position of the particle along the 
washboard corresponds to the phase difference across the 
junction, so that the speed of the particle as it moves 
in the washboard potential corresponds to the voltage 
across the junction. As it moves along the washboard, 
the particle is subject to a viscous damping force which 
is inversely proportional to the resistance shunting the 
junction. The strength of the damping can be character- 
ized by a quality factor parameter^^ Q = upRC, where 



R and C are the resist ance and capacitance shunting the 
junction and up = y^2eIc/hC is the angular frequency 
of smah oscillations at the bottom of the potential well 
at zero bias. Hysteretic IV characteristics are obtained 
for Q 3> 1 and phase diffusion obtained for Q ~ 1. 



1.2. Characteristic rates 
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where AC/g is the height of the energy barrier from 
a washboard potential minimum to the adjacent max- 
imum, a^ is a damping-dependent pre-factor and the 
quantities af^,uja, and AC/g are all current dependent, 
with iua, = ujp{l - [I/IcfY^^ and, close to /c, AC/g w 
^^-Ej(l — I/Ic)'^^^ where the Josephson energy Ej = 
hlc/2e. Combining Eqns. [l]and[2]gives a characteristic 
asymmetric distribution of switching currents for such 
junctions, as shown in Fig.[T^. 

In the underdamped regime, the mean switching cur- 
rent decreases as the temperature increases (Fig. [T]d) be- 
cause larger thermal fluctuations enable escape from the 
washboard minimum at a lower current. The width of the 
switching current distribution may be shown to depend 
on temperature as cr '--^ T^'^ . 

Thermal fluctuations can also cause retrapping of a 
particle which has escaped from a potential well. Ben- 
Jacob et al.^ obtained an analytic formula for the rate Fp) 
of this retrapping in the limit Q 3> 1. The retrapping rate 
is strongly dependent on the damping through Q, and is 
given by 
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where /r = hiQ) ~ AIq/ttQ. Rewriting this in the 
form Fj^ ^ exp (— At/j^/ZcT) defines an energy barrier 
AC/j^ for retrapping.^'* Eqn. [3] has been applied in the 
literaturei^ in the regime of moderate damping Q ^ 5, 
and we consider here in further detail application of the 
model in that regime. 

It is instructive to define a normalised current-ramp 
rate Fj = jd//dt. Eqn. [l]can then be rewritten 
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FIG. 1: (Color online) (a) Calculated underdamped thermally 
activated p{I) distributions at a number of temperatures. The 
parameters indicated are used for subsequent figures unless 
otherwise stated; these parameters might be typical for an 
IJJ.— Vertical lines (pink and blue online) show the boundaries 
of the distribution at 5 K, within which 99.99% of switching 
events occur. The thick vertical line (red online) shows the 
critical current, (b) Variation with temperature of the mean 
switching current (thick black line) and the top (blue online) 
and bottom (pink online) of the switching current distribu- 
tion. The lower and upper broken lines (red and green online) 
show respectively /(Fg = Fj) and /(Fj^ — lOFj). 



The rate of thermally activated escape from a mini- 



For small currents, Fg ^ Fj, so p{I) is small. As the 
current is increased towards the current /gj, at which 
Fg = Fj, the first quotient in Eqn. |3] increases and there- 
fore, as the current increases further, the numerator of 
the second quotient^^ begins to reduce from 1 to zero. 
The maximum in p{I) therefore occurs for Fg > Fj. 
The dashed lines in Fig. [T|d show the currents at which 
Fg = Fj and Fg = lOFj. The exact ratio Fg/Fj at 
the maximum in p(I) is temperature dependent: at 5 K, 
the peak in the switching current lies at a higher current 
than the current at which Fg = lOFj, whereas at 15 K, 
the peak in the switching current lies at a lower current 
than the current at which Fg = lOFj. 



THE MULTIPLE SWITCH-RETRAPPING 
REGIME 




FIG. 2: (Color online) Variation of characteristic rates with 
current. The different panels show the effect of variations in 
temperature and Q. 



As the current is increased from zero, the three char- 
acteristic rates Fg, Fj^ and Fj vary. Fig. [2^ shows the 
variation of these three rates when Q = 200, i.e., for an 
underdaniped junction. At very low currents, the retrap- 
ping rate is much larger than Fg and Fj. Also, since 
Fg ^ Fj, no escape events occur. When the current is 
increased to around /gj, an escape event becomes likely, 
but for / ^ /t7t the retrapping rate is very much smaller 
than the escape rate. Therefore retrapping is negligible 
in the case illustrated in Fig. [5^. As we will see, an im- 
portant current is the current /gp? at which Fg = Fjs . 
In Fig. Hi, /gj^ = 0.31 ^J^k and /gj^ < /gj. 

Fig. ^Bp shows the variation of the three characteristic 
rates for a more heavily damped junction. The escape 
rate is only weakly dependent on Q through the pre- 
factor a^ (Eqn. [J). The retrapping rate is exponentially 
dependent on Q (Eqn. [3]); it is much larger in Fig. ^Bp 
than in Fig. [5^ and Jgj^ > /gj. For currents / ^ /gj, 
the retrapping rate Fp; is now much larger than Fg . Es- 
cape events occur for / ^ /gj, but retrapping occurs 
shortly afterwards; the particle moves down the wash- 
board in fits and starts and the time-averaged voltage 
across the junction is non-zero — this state can be called 
a region of phase diffusion!^ As the current increases, the 
escape and retrapping rates become more and more sim- 
ilar, so there is a gradual increase in the time-averaged 
voltage. Fig. [3] shows, for the same values of Q and T 
as Fig. ^Bp, a simulation of jumps between the super- 
current (zero voltage) and running (resistive) states at 
three currents close to /jrjj^. In Fig. [3^, / < /gj^ and 
the junction spends most of the time in the zero-voltage 
state. At / '^ /gji (Fig. Wp)i escape events and retrap- 
ping events are expected in similar proportion and the 
junction spends a similar amount of time in the zero- 
voltage and escaped states. The time-averaged voltage 
across the junction becomes a significant fraction of the 
fully switched voltage, so an experiment is likely to mea- 
sure a switch event. As the current is increased further 
above igj^, any retrap event will be followed quickly by 
an escape event, so the junction spends almost all its 
time in the running state, as Fig.[3J: shows. For the junc- 
tion parameters corresponding to Fig. [2]3, the junction 
switches around /gpj > -^gj, so the switching current is 
greater than the switching current in the underdamped 
case. In other words, counter-intuitively, thermal fluctu- 
ations suppress the switching current less in the multiple 
switching-retrapping regime than in the conventional un- 
derdamped thermally-activated switching regimci^ 



3. TEMPERATURE DEPENDENCE, THE 
CROSSOVER REGIME AND T* 

Figs. ^Bp and c show the variation in the characteris- 
tic rates for two different temperatures with Q = 7. At 



R 
S 




< > 
0.5|xs 

Q=7 
T=20K 

(a) l<l,. 











3.59 



3.6 3.61 

current (|xA) 



3.62 



R 
S 


















< > 
0.5|^s 

(b) !==!,, 



















3.82 



3.83 
current (|xA) 



3.84 



R-i 



S "-I 



0.5)j,s 



(c) l>l 



ER 



3.94 



3.95 



3.96 
current (|a.A) 



3.97 



3.98 



FIG. 3; (Color online) Switching between zero- voltage super- 
current state (S) and the resistive running state (R) at three 
currents close to /go = 3.8417 /iA in a representative sim- 
ulation for the junction parameters shown in Fig. ^jp. The 
instantaneous voltage in the running state well above the re- 
trapping current is given by IR, where R is the relevant re- 
sistance shunting the junction. 



the lower temperature, 5 K, (Fig. [2j;) /gj > /gj^ so, for 
Fg '^ Fj, the retrapping rate is smaller than the escape 
rate. Therefore, retrapping after escape does not occur, 
and the conventional underdamped thermal activation 
behaviour is obtained. Conversely, at a higher tempera- 
ture, 20 K, (Fig. [IId), /gj < Iji]j{, and so there are mul- 
tiple escape and retrapping events, as described earlier. 
Note that /gp? is approximately unchanged as tempera- 
ture varies at constant Q. This is expected by inspection 
of Eqns. [5] and [31 Ignoring corrections of logarithmic 
order, AC/t7(/t7p) = AC/Ti(/g-D) determines Ivt>, where 
AC/g and AC/jj are both independent of temperature and 
hence -/gn is independent of temperature too. 




T(K) 

FIG. 4: (Color online) (a) Variation of standard deviation of 
switching distribution for Q = 7 with temperature, (b) Vari- 
ation of mean switching current with temperature for Q = 7. 
Shown are simulation results (open circles, with line to guide 
the eye, blue online), the full underdamped thermal (thick 
line, red online) and the underdamped thermal distribution 
truncated to above /g]^ (small closed circles, with black line 
to guide the eye). The dash-dotted line shows I^yi '^'^'^ ^* i'^ 
defined by the maximum in the simulated distribution width. 
(c) Variation with temperature of the mean number of es- 
cape events before a switch is counted in simulations. In the 
underdamped thermal case, a single escape event would lead 
to a switch. Above ri*- i , the distribution in this value is 
smaller than the symbols. Simulated distributions were based 
on 25000 switching events. 

The results of simulations of these dynamics over a 
broader range of temperatures are shown in Fig. 3] In 
Fig. 2^, the width of the distribution follows the con- 
ventional underdamped thermal behaviour {a ~ T"^'^) at 
lower temperatures, passes through a maximum and then 
falls at higher temperatures, matching experimental ob- 
servations. In the previous experimental reports ji^^iii 
a characteristic temperature T* was defined as the tem- 
perature at which the maximum in a{T) occurs. At low 
temperatures, a single escape event results in a switch 
being counted (Fig. Uj;) and the mean of the distribu- 
tion follows the conventional underdamped thermal be- 



haviour (Fig. [4)3). At around T*, the mean switching cur- 
rent flattens out and reaches an approximately constant 
value I « ^ER ^^-^^ above T* . For higher temperatures 
a significant number (^ 10'^ — 10'* above 25 K) of escape 
events occurs before a switch is counted. The shape of 
the switching distribution also changes as the tempera- 
ture is increased. Fig. [Sb-e shows that the shape departs 
from that shown in Fig. [TJ The skewness (the ratio of 
the third moment about the mean to the standard de- 
viation) gives a simple one-parameter description of the 
shape of the distribution; a symmetrical distribution has 
zero skewness. The skewness of the underdamped ther- 
mal distribution is around — 1 over the range of temper- 
ature shown in Fig. IH Fig. [5^ shows the variation of 
the skewness of the simulated distribution around T* . 
The skewness of the distribution begins to depart from 
its thermal value somewhat below T* , becoming progres- 
sively less negatively skewed and then positively skewed, 
passing through a maximum and then beginning to level 
out at around the same temperature as the width begins 
to level out and as the mean switching current levels off. 

From these simulations, we identify three different 
regimes of behaviour. At low temperatures, conventional 
thermal underdamped behaviour is observed. At some 
higher temperature below T*, the skewness of the dis- 
tribution, and in detail also the width and the mean, 
depart from the underdamped thermal values. Above 
this temperature, the skewness and width vary rapidly. 
At a higher temperature, there is a crossover to a dif- 
ferent regime in which the mean switching current is ap- 
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FIG. 5: (Color online) (a) The skewness of the switching dis- 
tribution, negative for conventional underdamped distribu- 
tions, becomes positive around 2T- v . Points and lines as 

described for Fig. |3] (b)-(e) Switching distributions at se- 
lected temperatures. Note the scales on these insets are the 
same, with the scale on the horizontal current axis running 
between 3.5 /xA and 5 /xA and the scale on the vertical p{I) 
axis running from to 10 (/lA)"'^. Simulated distributions 
were based on 25000 switching events. 



proximately constant and the skewness and width are 
slowly decreasing as the temperature is increased. To 
describe this behaviour, we label the two boundaries 
between these three regimes T-,* and Tv*. , , where 

low high' 

To arrive at a quantitative definition for T,* , we note 
that retrapping only has a significant effect on the dy- 
namics when there are escapes at currents / < /jT;pj. 
Therefore, for parameters where there are no escapes 
for / < /jr;p>, the switching distribution does not de- 
part from the conventional underdamped thermal distri- 
bution; we define Ti* quantitatively as the temperature 
at which /jt;p) coincides with the bottom of the conven- 
tional thermally activated underdamped switching dis- 
tribution (see also the lower line (pink online) in Fig.[T]), 
where we define the bottom I-^ and top I\^ of the distri- 
bution by /o^p(/)d/ = /p and /5*^p(/)d/ = /p, where 

< /p <C 1, with /p = 0.0005. Fig. [5^ shows a simulated 
switching-current distribution at TT and a comparison 
of Fg, Fj} and Fj as a function of current. 

As the temperature increases, the current /gj de- 
creases. For temperatures T > T-,* , escape events occur 
for /gj ^ I ^ ^ER ^"^^ ^^^ followed by retrapping events; 
they do not result in the count of a switch. Fig.[5jD shows, 
for a temperature T > T^ , a comparison of the charac- 
teristic rates and also shows that the simulated switching 
distribution begins to depart from the conventional un- 
derdamped thermally activated switching distribution. 

For T > Ti* , escape events leading to switching only 
occur at currents / > ^eR- When a part of the un- 
derdamped thermal distribution lies at / > ^ER' '^^^ 
might naively expect that the width, mean and shape of 
the switching distribution would be approximately the 
same as the width, mean and shape of the part of the 
underdamped thermal distribution lying at / > ^eR- 
The mean, width and skewness of the part of the un- 
derdamped thermal distribution lying above -^eR ("the 
truncated distribution") are shown in Fig. S] and Fig. [Sj 
The mean in the simulations closely matches the mean 
of the truncated underdamped thermal distribution, and 
the temperature dependence of the width and the skew- 
ness in the simulations also follow the respective vari- 
ations in the width of the truncated distribution, up to 
temperatures approaching the temperature at which ^eR 
coincides with the top of the conventional thermally acti- 
vated underdamped switching distribution (see also the 
upper line (blue online) in Fig. [T|). We define this lat- 
ter temperature as Ti*. , (see Fig. [6};). The shape of 

the part of the underdamped thermal distribution with 
/ > -^ER largely determines the shape of the switching 
distribution for T,* < T < Tj*. 1 ,i.e., the conventional 

thermal behaviour is only followed below T,* , but the 
variations in the mean and width of the distributions re- 
main analytically determinable for temperatures up to 
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FIG. 6: (Color online) Variation of characteristic rates with 
current at three temperatures for Q = 7, along with the 
switching distribution. The heavy (blue online) points show 
the simulated switching distribution. The lower black curve 
shows the underdamped thermally activated switching distri- 
bution. The broken black line depicts the current /]7;r . (a) 
T — ri* — 9.8 K; /]t;r coincides with the bottom of the ther- 
mal distribution, (b) T = 13.0 K. (c) T = T^^„^ = 16.0 K; 
-^ER coincides with the top of the thermal distribution. Sim- 
ulated distributions were based on 100000 switching events. 



The departure of /gw from the conventional thermal 
activated behaviour above T,* and the approach to an 
asymptotic value is in agreement with experimental re- 
sults reported by Franz et ai, in which a plateau in 
^sw(T') is observed above a crossover temperature, fol- 
lowed by a fall at higher temperatures. In one sam- 
ple, however, they observe an increase in /gw above the 
crossover temperature. These observations may be ex- 
plained by considering temperature variation in /c and 
Q. The Ambegaokar-Baratoff relation^ suggests a sig- 
nificant decrease in /c for T > Tc/2. A reduction in /c 
as T increases would lead to a fall in the switching cur- 
rent, whereas a decrease in Q could lead to an increase 
in the switching current. One might expect a reduction 
in Q at higher temperatures as quasiparticle conductiv- 
ity increases. This change will be sample dependent; for 
a typical IJJ sample, the resistance may fall by a third 
from low temperature to Tq- However, if the shunt resis- 
tance is dominated by the environmental impedance, Q 
is expected to be approximately T-independent. 

Since, as the temperature increases above T,* , the 
simulated distribution becomes progressively less neg- 
atively skewed and becomes positively skewed around 
T^*. 1 , the temperature variation of the skewness pro- 
vides a straightforward experimental way to determine 
r,* and Ti*. , . As far as we are aware, no systematic 

experiments investigating variations in the shape of the 
switching distribution around T* have yet been reported. 

For T > TZ- u, the behaviour is not associated with 
^ high' 

the underdamped thermal distribution; we will not dis- 
cuss this behaviour in detail. Fig. HJd shows that, for 
T > Tj*. 1 , the mean switching current approaches /gj?. 

Since /gj? is a function of Iq and Q, determination of 
/gp> from the switching current at T > Tj*. , provides 

an additional experimental probe for the determination 
of /c and Q, and a consistency check for derivation of Q 
from either T* or from a ratio of the switching and return 
currents. In addition, the expectation that the mean of 
the switching distribution at temperatures above TV*. , 

does not vary with temperature for constant Q means 
that this measurement might be used as an experimental 
probe for whether Q is varying with temperature. We 
return later to discuss how the presence of frequency- 
dependent damping affects this plateau. 

For T > T^ , the behaviour shows some similarities 
to the extensively studied phenomenon of phase diffu- 
sion, but also some differences. The time-averaged volt- 
age across the junction is finite below the switching cur- 
rent, as expected for conventional phase diffusion. How- 
ever, for / < /gj, the time-averaged voltage across the 
junction is identically zero with a high probability.— 
This behaviour contrasts with a phase-diffusion regime 
considered by Ivanchenko and Zil'bermanl" and others 
in which a finite phase-diffusion resistance persists even 
at currents / ^ 0. For / > /gj, the time-averaged 
junction voltage is finite, but this phase-diffusion regime 
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FIG. 7: (Color online) Calculated variation in Ti™ and 2T- n 

with ramp rate, with /p=0.0005 and T/^^, the temperature 
at which Fg, Fr and Fj are equal at some value of current. 
Also shown is the temperature where the maximum width of 
the truncated (/ > Ij^Ji) distribution lies. 



also differs physically from that analysed by Ivanchenko 



and Zil'bermar 



Ivanchenko and Zil'bernian consid- 



ered that each escape vifould lead to a phase shift of 
only 27r, whereas in our model the phase shifts are of 
order ujp/T-n S> 1 — Kautz and Martinis previously 
showed that the presence of multiple-27r phase shift- 
escape events has an appreciable effect on the phase- 
diffusion voltagCiii This is associated with the time ~ 
l/wp needed after the energy barrier is exceeded for the 
instantaneous voltage to increase from zero to its steady 
value. In previous experiments on moderately damped 
junctions, Krasnov et al^ did not observe any phase- 
diffusion voltage above T* until well above T*. This ap- 
pears to conflict with our understanding from our anal- 
ysis. However, the explanation might simply be that the 
phase-diffusion voltage was too small to be measurable. 
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FIG. 8: (Color online) (a) Variation of width with tempera- 
ture around T* for various ramp rates, (b) Variation of width 
with ramp rate for various temperatures close to T* . Each 
datapoint corresponds to a simulated distribution based on 
1000 switching events. Departures from a smooth trend are 
visible as a result of statistical fluctuations. Lines are guides 
to the eye. 



3.1. Ramp-rate dependence 

Since the crossover temperatures are dependent on the 
shape of the thermal distribution, and the shape of the 
thermal distribution is dependent on the current-ramp 
rate, the crossover temperatures are also dependent on 
the current-ramp rate. Figure [7] shows the variation with 
current-ramp rate of calculated values of T,* and 2T- i 

and the temperature T* at which the maximum in the 
width occurs. All these values were determined from the 
truncated thermal distribution. Also shown is T^r^r. , 
the temperature at which /gj? = /gj, in this analysis 
the most natural definition for a single crossover tem- 
perature. For a ramp rate d//di = 10~^ As~^, we find 



T, 



low 



8.2 K and T^- . 
high 



12.2 K, whereas for a ramp 



rate d//dt = lO^^ As"\ we find T^^^ = 10.8 K and 
T|*. 1 = 18.9 K. Thus, as the simulations reported in 

Fig. [8] show, varying the ramp-rate can have a significant 
effect on T*; this dependence was not recognised in pre- 
vious reports of the crossover. Fig. [5^ shows simulations 
of the variation in the width of the switching distribution 
with temperature for a number of current-ramp rates. In 
the range shown, for dl/dt = 10^"* As^^ the temperature 
is above T,* , whereas for dl/dt = 10~^As~"'^ the tem- 
perature ranges from below T,* to above T,*. , . This 

difference between the two ramp rates leads to a marked 
difference in the temperature variation of the width. 

Experimentally, probably the most straightforward 
measurement to make to investigate these effects would 
be to keep T fixed and vary d//dt. Fig. [5)3 shows a 
simulation of this procedure — there may be a marked 



difference in the variation of the width with ramp rate 
depending on the temperature of the measurement. For 
T = 9 K, the temperature is at or below T,* for the 
whole range of ramp rates, whereas for T = 18 K, the 
temperature is at or above Tj*. i for the whole range of 

ramp rates. For intermediate temperatures, increasing 
the ramp rate from lO""* As~^ to 10~^ As~^ moves the 
T* values so that the temperature is close to T-i*. , at the 

lowest ramp rate and close to T,* at the highest ramp 
rate. This ramp-rate dependence of the width could be 
used to determine Q. As the ramp-rate is varied at a con- 
stant temperature, Q remains fixed, but T* varies and so 
the width of the distribution varies, particularly when 
T* becomes close to the experimental temperature. This 
allows Q to be determined, at each experimental temper- 
ature, by fitting to simulations. 



4. THE RETURN CURRENT h, HYSTERESIS 
AND FREQUENCY-DEPENDENT DAMPING 

In previous publications^ii^ it was noted that, although 
the temperature at which the hysteresis in the IV char- 
acteristic disappeared was around T*, there was some 
difference between the two values. Motivated by this dis- 
crepancy, we consider here in more detail the variation 
of hysteresis around T* . 

In the RCSJ model in the absence of fluctuations, as 
an applied current is ramped down from Iq towards zero, 
return from the quasiparticle branch occurs at a current 
Ij^ « 4:1c /ttQ for Q ^ 3.1- Thermal fluctuations lead to 
retrapping when Fj^ ~ Fj, at currents /r > /j^. As Fig. [9] 
shows, there is some distribution in the value at which 
return occurs, and the mean and peak of the distribu- 



tion lie below the current at which F 



R 



Fj. The width 



of the return distribution may be shown to vary with T 
as cr ~ ri/2^ Experimentally, measurements of return 
distributions are more likely than switching-distribution 
measurements to be affected by heating and this com- 
plicates analysis of the temperature dependence. To our 
knowledge, the only report in the literature of an exper- 
iment in which the distribution of return currents was 
measured has been given by Castellano et al^ 

A measure of the hysteresis of the junction is given by 
Ir/Isw and, experimentally^*, Ir/Isw is sometimes used 
to infer Q through the approximation /r/^sw ~ -^r/-^c — 
-i/TrQ. However, for moderately damped Josephson junc- 
tions, the ratio Ej/kT is often sufhciently small that 
thermal fluctuations cause significant departures of /j- 
from /j5 and Isw from Ic- Fig. [10] shows the varia- 
tion of Ir/Isw with area at selected temperatures, as- 
suming junctions remain in the underdamped regime. 
The junction parameters chosen might be appropriate 
for intrinsic Josephson junctions. For large-area junc- 
tions. It/ Isw ~ ^r/^c- For junctions with area ~ 1 ^m^, 
Ii/Isw is significantly larger than I^i/Ic ^nd so in this 
case identifying Ir/Isw with Iji/Ic to infer Q will give 



a significant underestimate of Q. The inset in Fig. [TOl 
shows the difference between the crudely inferred Q and 
the true Q at 4.2 K: the crudely inferred Q is less than 
the actual Q. If measurements at higher temperatures 
were used, the discrepancy would be larger. Therefore 
it is important to account for the reduction by thermal 
fluctuations of the switching and return currents. We em- 
phasize that, since we are treating Q as a temperature- 
independent quantity, this reduction in the hysteresis is 
purely a thermal effect. For junctions with area 0.1 /xm^, 
it can be seen by comparing Figs. [T] and [5] that the un- 
derdamped return current and underdamped switching 
current become similar around 14 K, and at higher tem- 
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FIG. 9: (Color online) (a) Thermally activated return distri- 
butions at various temperatures, for Q = 7. The distribution 
was constructed from the retrapping rate in Eqn.[3]in the same 
way that Eqn. [2]was used to construct the p{I) distribution 
shown in Fig. [T] Vertical lines (pink and blue online) show 
the boundaries of the 5 K distribution, within which 99.99% 
of switching events occur. The thick vertical line (red online) 
shows the return current in the absence of fluctuations, (b) 
Variation with temperature of the mean return current (thick 
black line) and the top (pink online) and bottom (blue online) 
of the return current distribution. 
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FIG. 10: (Color online) Area dependence of the ratio of mean 
return current and mean switching current for a quasi-dc mea- 
surement d//dt = 10~^ As~^. The dashed line shows the 
fluctuation-free value Ir/Isw = -^r/^c = 4:/ttQ. Inset: Varia- 
tion of inferred Q — AIsw /ttIt at 4.2 K with Q for two junction 
areas. The dashed line shows Qi^ferred ~ '3 actual- 




FIG. 11: (Color online) Temperature dependence of the ra- 
tio of mean return current and mean switching current for 
d//df = lO"'^ As~^. The gray line (red online) shows the 
variation derived from the conventional thermal activation 
(Figs, [l] and |9]) for (5 = 7. Black points show the dependence 
obtained from simulations including both escape and retrap- 
ping. The dashed horizontal line shows the value at which 
hysteresis disappears Ir/Isw = 1. Inset: Simulated switch- 
ing distribution with current ramping up (heavy black curve) 
and simulated return distribution with current ramping down 
(heavy gray curve — red online), for T=14.5 K, demonstrat- 
ing hysteresis in the switching and return currents. For com- 
parison, the corresponding conventional thermally activated 
distributions are also shown (black and gray (red online) thin 
lines) . The broken vertical line shows the position of /jt;r and 
the dotted line (blue online) shows (the unrescaled) P^Pjifiil) 
from Ref. JA for comparison. Compare with the fluctuation- 
free values, 7c = 10/xA and Jj^ = AIc/nQ = 1.82 /lA. 



peratures the underdamped return current exceeds the 
underdamped switching current. This would imply that 
there is a current range Isw < I < It where neither 
the zero-voltage branch nor the resistive branch is sta- 
ble — this is indeed a feature of the behaviour we are 
describing in the temperature range T > Ti* (for ex- 
ample, see Fig. [3]). Above T,* , as we have discussed, it 
is necessary to consider, in addition, the effects of mul- 
tiple escape and retrapping. Fig. [TT] shows the variation 
with temperature in the ratio Ir/Isw of these mean val- 
ues. This ratio was obtained from simulations including 
the effects of both escape and retrapping as the current 
is ramped in cither direction. At low temperatures we 
find Ir/Isw ~* '^/iiQ. As the temperature is increased, 
even well below T* , Ir/Isw departs significantly from 
i/nQ. Around T* , the distribution departs from its con- 
ventional thermally activated behaviour. Hysteresis is 
still present for T,* < T < T^. , , but to a decreasing 

extent as T is increased. The inset of Fig. [TT] shows the 

switching and return distributions for a temperature 14.5 

K where T,* <T < T,*- ,. The thick black curve shows 
low high 

that the mean and peak of the switching distribution lie 
above It?t} ■ The thick gray curve (red online) shows that 
the mean and peak of the return distribution lie below 
-^ER (^^^ width is related to the shape of the conventional 
return distribution). A difference in the mean switching 
and return currents, and so some hysteresis, persists even 
though the underdamped return current exceeds the un- 
derdamped switching current, which might be thought to 
imply the absence of hysteresis. The persistence in hys- 
teresis above T* in the experiments of Refs. [ij and [^ 
is likely to be attributable to the distinction between 
T* and Tu- , — we expect hysteresis to persist up to 

T 



For T > Ty^. , , both escape and re- 

1 and 



and Tr. , — 
high 

'" high 
trapping occur close to /gR, so that Ir/Isw 
hysteresis in the IV characteristic is small. 

We would like to note that the probability of a switch 
being counted is not the same as the probability of a 
single escape not followed by a retrapping event, since 
a switch may be preceded by many escape-retrapping 
events. This difference was not appreciated in the quan- 
titative analysis in Ref. Il4i in which the latter quantity 
— although much smaller than the total probability of 
a switch — was evaluated as a function of the current 
of the initial escape (compare the dotted curve and the 
thick black line in the inset of Fig. fTTj) and then rescaled— 
and fitted to experimental data. 



4.1. Frequency-dependent damping 

In our treatment so far, we have been assuming that 
the damping Q is frequency independent. Different re- 
trapping and return behaviour may arise when the damp- 
ing is frequency dependent and we now turn to dis- 
cuss these differences. If the damping of the system is 
frequency-dependent — as is likely to be the case un- 
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less isolation resistors close to the junction are included 
or the shunt resistance is much less than the free-space 
impedance — the system will be characterized by much 
stronger damping shortly after escape than in the steady 
running state. This means that overdamped behaviour 
at escape might lead to multiple escape-retrapping be- 
haviour shortly after an initial escape but, once the junc- 
tion has been in the running state for some time, it be- 
comes much less strongly damped and might be charac- 
terized by the conventional underdamped dynamics and 
so unlikely to be retrapped. 

In their simulations, Kautz and Martinisii consider 
a crossover to low-frequency damping once the junction 
has been mostly in the running state over a timescale 
l/i^C7 where vc is a characteristic frequency, a factor of 
10^ — 10^ smaller than the plasma frequency. In their 
simple model, which captures the qualitative features of 
the frequency-dependent-damping behaviour, well above 
Vc the damping is Qi and well below uc the damping is 
Qo < Qi- 

In the simulations presented here so far, we counted a 
switching event if the junction was mostly in the running 
state over a certain time period. Experimentally, that 
time period might be determined by the response time of 
the measurement electronics. In the case of frequency- 
dependent damping, the time period is set instead by 
l/j^c, since once the junction has been mostly in the run- 
ning state for the characteristic time l/i'c and becomes 
characterized by the much lighter damping Qi, Fj^ falls 
and retrapping becomes unlikely. This time is likely to be 
much shorter than the response time of the electronics. 

To model, through simulations, the effect of a decrease 
in the damping at low frequencies, we neglect retrap- 
ping once the junction has been in the running state 
for a time I/vq. Fig. [T2] shows the variation in the 
mean and width of the switching distribution around 
T* as the characteristic time I/vq varies. Note that 

the crossover temperatures T,* and T,*. , are essen- 
'■ low high 

tially independent of the measurement time period, since 
they are set by the crossover between the extremes of 
the conventional underdamped distribution and by /gj^, 
all of which are essentially independent of the measure- 
ment time period. Fig. [T^] shows that a decrease in the 
characteristic timescale leads to a decrease in the mean 
switching current and to an increase in the width above 
T* . The width and mean of the switching distribution 
for T > 2T- .1 therefore provide an indication of the 

crossover frequency vc- 

This variation with vc rnay be understood by consid- 
ering the relative sizes of Fg and Fp> , in comparison to 
Vc, around the current at which switching occurs. In 
the multiple escape-retrap regime, for small vc such that 
Vc <S Fg, Fp> , the fraction of time in the running state /r 
is to a very good approximation /r,eq = ^^/(Fg -I- Fj^) 
and so switching is likely for Fg = Fp> , i.e., for I = /gpj . 
(For the simulations presented in Section [31 the same 
approximation holds for the temperatures of interest 
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FIG. 12: (Color online) Variation of distribution characteris- 
tics with the timescale l/j^c for Q — 7. Main panel; Mean 
switching current; the dash-dotted line shows /g]^. Inset; 
Switching distribution width. Lines joining points are guides 
to the eye. The lines without points (red online) show the 
conventional underdamped thermal behaviour. Each data- 
point corresponds to a simulated distribution based on 1000 
switching events. 



T « T*, since the time-period t = 10"^ s > l/Fg, l/Fj^ 
— compare Fig. [HI) However for larger vc such that 
1/i/c < (1/Fg-|-1/Fj|), large fluctuations of /r away from 
/r,eq occur. These fluctuations make switching likely at 
/ < Iji^Yi ^^'^ ^^^° increase the width of the switching 
distribution. The larger vc is, the smaller Igw is likely to 
be. For example, in Fig. [5^, for l/i'c = 10~^ s, the es- 
cape to the running state at around 3.603 fiA would last 
long enough to cause a switch, although /gj^ = 3.8417 
fiA. 

The presence of frequency-dependent damping also has 
a marked effect on the hysteresis of IV measurements. 
For an initial escape as the current is ramped up from 0, 
the system is characterized by the high-frequency damp- 
ing Qo- In contrast, when the current is ramped down 
from /c I the system is characterized by the low- frequency 
damping Qi. The more underdamped behaviour for re- 
turn means that hysteresis persists to much higher tem- 
peratures than for a system with frequency-independent 
damping Qo- 



4.1.1. Application to previous work 

We now briefly discuss previous work in relation to 
frequency-dependent damping. The work of Krasnov et 
a/.^— lis, is largely on samples with resistances ~ 1 fi which 
are therefore likely to be characterized by approximately 
frequency-independent damping, although the larger re- 
sistance of the IJ J samples suggests they may be affected 
by frequency-dependent damping, so that the damping 
characterizing retrapping soon after escape is not the 
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same as the damping extracted from the hysteresis in 
the IV characteristics. 

In their paper, Kivioja et at— were considering mostly 
the switching dynamics, which are determined by the 
high-frequency damping, although the retrapping current 
/m is determined by the low-frequency damping. They 
did not include frequency-dependent damping in their 
modelling and they used the high-frequency i? as a fit- 
ting parameter. They expect that the maximum in the 
width of the distribution to occur at T^, when r-g(/ni) 
corresponds to their experimental timescale, essentially 
equivalent to our Fj. This contrasts with our expec- 
tation for frequency-independent damping that, at T*, 
Isw = /eR- Si'^ce T-^ilm) < ^^{Ier), the crossover 
temperature T* < Tj. Kivioja et al. extracted a value 
for Q from their experimental results which is therefore 
larger than would be extracted if thermally activated re- 
trapping had been included. This may explain the dis- 
crepancy between their extracted Q = 4.4 and the (3 = 4 
suggested by the nominal values of their experimental 
parameters. 

Mannik et al}'' obtained values for the probability of 
retrapping from numerical Monte Carlo simulations and 
included frequency-dependent damping. This model has 
the advantage that it is able to account for the presum- 
ably initially increased rate of retrapping as, in the wash- 
board analog, the particle first accelerates after escape. 
However, the results are less straightforward to analyse. 
The authors expressed the net escape rate as a sum of 
the probabilities of multiple escape-retrap events, related 
to the thermal escape rate (our Eqn. [J) and the calcu- 
lated retrapping probability. Although in detail the as- 
sessment of individual escape probabilities, a nontrivial 
problem, is oversimplified and relies on strictly inconsis- 
tent approximations^ the model captures at least qual- 
itatively well the effect of retrapping on the switching 
current statistics. Mannik et al. found good agreement 
between their model for the net escape rate and their ex- 
perimental results. Their treatment of the probability of 
retrapping after escape as a time-independent quantity 
contrasts with the model of Ben-Jacob et al^, in which 
retrapping is modelled by a rate (Eqn. [3]) and so with 
a probability increasing linearly with time spent in the 
running state. The treatment of this probability as time- 
independent by Mannik et alJ^ is successful because, in 
the model of frequency-dependent damping which they 
used, the time over which retrapping can occur is « l/j/c. 
Once the particle has been in the running state for a time 
> l/i^c ~ 1/<j-'Pj low- frequency (under)damping applies 
and retrapping is unlikely. This expectation is borne out 
by their simulations, in which they observed retrapping 
events only times < lOO/wp after escapeJ^ 



damped Josephson junctions. We emphasized that there 
is a regime in which the junction repeatedly escapes to 
and retraps from the running state and demonstrated 
through the use of simulations that for some choices of 
parameters, the number of escapes and retraps during a 
single current ramp to an eventual switch into the run- 
ning state may be very large (~10000). The multiple 
escape-retrapping regime, with a large number of escapes 
of duration ~ l/fj^, is intermediate between the under- 
damped regime in which a single escape leads to switch- 
ing, and the overdamped phase-diffusion regime in which 
a very large number of escapes of very short duration 
~ 1/^p may occur. 

By examining the region around the crossover in the 

temperature dependence of the width in more detail, we 

showed that the crossover is, in detail, described by not 

one but iiwo crossover temperatures T,* andT,*. , . The 

^ low high 

variations in the mean and width of the switching dis- 
tribution (in the intermediate regime between the two 
crossovers) are largely determined by the shape of the 
thermally activated switching distribution and this shape 
therefore also determines the temperature of the maxi- 
mum in the width, the quantity usually identified as the 
single crossover temperature. We showed that the shape 
of the switching distribution, parametrized by the skew- 
ness, indicates T-,* and Ti*. , . We introduced a perti- 
nent rate Fj for understanding the dynamics; we showed 
that the details of the frequency dependence of the junc- 
tion damping should affect measured values of the mean 
and width of switching distribution and weakly affect the 
crossover temperatures. 

We showed that the characteristic temperatures T*, 
T,* and Ti*. , are all dependent on the current ramp- 
rate and therefore are not uniquely determined by mea- 
surements at a particular ramp rate, and in addition that 
there is some dependence of the behaviour around and 
above T* on any frequency dependence of the damping. 

We also considered the process of return to the super- 
current state as the current is ramped down in the pres- 
ence of thermally activated retrapping events and the im- 
plications for measurements of hysteresis in moderately 
damped Josephson junctions. We found that some hys- 
teresis is expected to persist above T* , to Tj*. , , even 

in junctions with frequency-independent damping. This 
suggests a resolution of the issue of hysteresis somewhat 
above T* in previous reports. 
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